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Particle-in-Cell simulations of collisionless magnetic reconnection with a guide field reveal for the 
first time the three dimensional features of the low density regions along the magnetic reconnection 
separatrices, the so-called cavities. It is found that structures with further lower density develop 
within the cavities. Because their appearance is similar to the rib shape, these formations are 
here called low density ribs. Their location remains approximately fixed in time and their density 
progressively decreases, as electron currents along the cavities evacuate them. They develop along 
the magnetic field lines and are supported by a strong perpendicular electric field that oscillates 
in space. In addition, bipolar parallel electric field structures form as isolated spheres between the 
cavities and the outflow plasma, along the direction of the low density ribs and of magnetic field 
lines. 



INTRODUCTION 

Magnetic reconnection is the topology reorganization 
of disconnected magnetic regions into a new configura- 
tion with concurrent conversion of magnetic energy into 
plasma kinetic energy. The magnetic energy is released in 
form of plasma heating and acceleration: plasma (inflow 
plasma) drifts toward the region where magnetic recon- 
nection takes place and expelled from it forming the re- 
connection jets (outflow plasma). The surfaces between 
the new magnetic regions, dividing the inflow and out- 
flow plasmas are called separatrices. Separatrices are thin 
boundary layers between plasmas with different proper- 
ties (density, temperature, and flow pattern) [T]-[3] and 
subject to instabilities, wave activity and turbulence [4]- 
[6]. One of the distinctive features of collisionless mag- 
netic reconnection is the formation of localized low den- 
sity regions along the separatrices [71 [8] . These depletion 
areas are called cavities. 

One of the most common and simple reconnection con- 
figurations is that of a sheared magnetic field layer, with a 
superimposed uniform magnetic field (called guide field) 
perpendicular to the plane of the reversing magnetic 
field (the reconnection plane). The configuration pro- 
vides a simple setup, but still realistic enough to compare 
simulation results with satellite observations of magne- 
topause [9UTT] and magnetotail [12]. In addition, this 
configuration results a useful approximation of magnetic 
field lines in nuclear fusion devices, such as the case of 



the tokamaks and reversed field pinch machines, where a 
toroidal field is present. In presence of guide field, sepa- 
ratrices display anti-symmetric features with the respect 
to the X point: the cavities develop only along two sep- 
aratrices, while regions of increased density form along 
the other two separatrices [6l[l3l[l4]. In addition, strong 
electron beams move toward the X line along the cav- 
ities separatrices. Streaming instabilities, such as the 
Buneman, two-stream, and lower hybrid instabilities are 
triggered by electron beams in the cavity regions. As 
macroscopic signature of these instabilities, bipolar elec- 
tric field structures are detected in simulations USHII], 
magnetospheric observations ^8| 119J and laboratory ex- 
periments [20 . 

The cavities were first identified in simulations by Shay 
et al. [21] , as one of the characteristic signatures of Hall 
magnetic reconnection. The decoupling of ions and elec- 
trons leads to an electron current system surrounding 
the magnetic separatrices, as suggested by the theory of 
Hall reconnection [22l |23] . The net parallel electric field 
strongly accelerates the electrons in this region, thus re- 
ducing the local density [24] . 

Previous studies focused on analyzing the structure 
and scaling properties of the cavities in two dimensional 
Particle-in-Cell simulations. It was found that cavities 
size scales as the ion skin depth, and their density is in- 
sensitive to the simulation mass ratio [14]. In addition 
the simulation results were confirmed by direct obser- 
vations of cavities in the magnetosphere by the Clus- 



ter spacecrafts [71 |8]. The majority of previous stud- 
ies mainly regarded to the two dimensional structures 
of the cavities, while the three dimensional structure 
was not studied in detail. This work focuses on the 
study of magnetic reconnection separatrices and three- 
dimensional cavity structure, presenting two main re- 
sults. First, it reveals the three-dimensional structures 
of the cavities. They appear as low density layers. Re- 
gions, characterized by further lower density, here called 
low density ribs form within the cavities. They develop 
along the magnetic field lines and are supported by strong 
perpendicular electric fields. Second, it shows that bipo- 
lar electric field signatures are present as spherical for- 
mations. They are aligned with the low density ribs and 
magnetic field lines, and they are located between the 
cavity and the outfiow plasma. 

This paper is organized as follows. The first section 
describes the simulation set-up and the parameters in 
use. The second section presents the evolution of three 
dimensional guide field magnetic reconnection, the cav- 
ity and low density ribs formation, and the presence of 
bipolar parallel electric field structures in their proximity. 
The third section discusses the results and the possible 
origins of three-dimensional structure of the cavities. Fi- 
nally, the conclusion section summarizes the results. 



A uniform guide field Bg = Bq along the z direction is 
added to the magnetic field profile and a perturbation of 
the z component of the vector potential, SAz^ is applied 
to initiate magnetic reconnection: 

dAz = A^o cos{27rx/LA)cos{7ry/LA)e-^''^^y^^^^^ , (3) 

with La = 10 A. 

The particles are initialized with a Maxwellian veloc- 
ity distribution. The electron thermal velocity Vthe/c = 
(Te/(mec2))^/^ = 0.045, while the ion temperature T^ = 
5 Tg. The simulation time step is UpiAt = 0.125. In 
this simulation set-up, ujpi/Vtci = c/Va = 103, where 
Qci = eBo/rrii is the ion cyclofrequency and Va is the 
Alfven velocity calculated with B = Bq and n = uq. The 
ion Larmor gyro-radius is 0.65 di. 

The grid is composed of 256 x 192 x 128 cells, resulting 
in a grid spacing Ax = Ay = Az = 0.08 d^. In total 3 x 
10^ computational particles are in use. The boundaries 
are periodic in the x and z directions, while are perfect 
conductor and refiecting boundary conditions for fields 
and particles respectively in the y direction. Simulations 
are carried out with the parallel implicit Particle-in-Cell 
iPICSD code i26j. 



RESULTS 



SIMULATION SET-UP 

Simulations are carried out in a three-dimensional sys- 
tem, where an Harris current sheet equilibrium is ini- 
tially imposed [25|. The z coordinate is taken along the 
Harris sheet current, while the x-y plane is the reconnec- 
tion plane. The x, y and z simulation coordinates corre- 
spond to the —xqsmi ^gsmi and —ycsM coordinates in 
the Geocentric Solar Magnetospheric (GSM) system and 
to Earth-Sun, North-South, and dawn-dusk directions in 
Earth magnetotail. 

The plasma density profile is initialized as: 



n{y) = no cosh 



y - Ly/2 

A 



Ub 



(1) 



The peak density no is the reference density, while the 
background density n^ = 0.1 no; L^ x Ly x Lz = 20 di x 
15 di X 10 di are the simulation box lengths and A = 
0.5 di is the half- width of the current sheet. The ion 
inertial length is di = c/ujpi^ with c the speed of light in 
vacuum, ujpi = y^^Tmoe^Jnii the ion plasma frequency, e 
the elementary charge, and nii the ion mass. The electron 
mass is mg = m^/256. 

A magnetic field in the x direction is initialized to sat- 
isfy the Harris equilibrium force balance: 



Bx{y) = ^otanh 



y - Ly/2 



(2) 



The evolution of magnetic reconnection is well depicted 
by the history of the electron current sheet, shown in a 
volume plot at different times in Figure 1. Different col- 
ors represent the intensity of electron current density: 
the green and blue colors identify high and low current 
density areas respectively. The initial thin current sheet 
breaks in the region where the magnetic reconnection was 
triggered by the initial perturbation, and separates in two 
parts moving out mainly along the x direction. It is im- 
portant to note that the presence of the guide filed alters 
the direction of the reconnection jets: their directions are 
not purely in the x direction, as in the case of antiparallel 
reconnection, but they are tilted by the guide field [12]. 
As result of the reconnection dynamics, thin current sur- 
faces, the separatrices, develop between the infiow and 
outfiow plasmas around time Qcit = 12.1. Finally, the 
separatrices become unstable at later times with the ap- 
pearance of filamentary structures along them. 

A new magnetic field topology is created from two dis- 
tinct initial magnetic regions as a result of reconnection. 
Figure 2 shows first a snapshot of the magnetic field lines 
at time Qcit = 12.1 in panel a, and then an isosurface 
plot of the reconnection electric field Ez from two differ- 
ent points of view in panels b and c. The reconnection 
electric field Ez reaches the peak value of 0.35 BqVa/c in 
a cylindrical region around the X line. The value of the 
reconnection electric field is in good agreement with the 
value obtained in two dimensional simulations of kinetic 
guide field reconnection [27H29] . 




FIG. 1. Volume plot of the electron current density at dif- 
ferent times. The initial Harris current sheet breaks in two 
reconnect ion jets, moving away from the x line. The electron 
current develops filamentary structures along the separatrices 
at later times. 
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As magnetic reconnection evolves, areas of decreased 
density regions, the cavities, develop along separatrices. 
Figure 3 shows a plot of the electron density evolution 
in an area around the X line {x = 4.8 di — 14.5 di^y = 
4 di — 11.2 di) on the reconnection plane at z = Lz/2. 
The electron density is normalized to the peak current 
sheet density no, and the maximum value of the color- 
bar is fixed to 0.3 no to investigate the low density re- 
gions. The typical density patter found in other two di- 
mensional Particle-in-Cell simulations is recovered. The 
density configuration is asymmetric with respect to the 
X point. The top left and bottom right quadrants of the 
four panels show the cavities, the low density regions in 
blue color in the contourplot. The cavities form at time 
Vtcit = 12.1, they are perturbed by a wave structure at 
time Qcit = 13.3 and are strongly distorted after time 
flcit = 14.5. Small spatial scale low density ripples de- 
velop along the cavities. These structures are similar in 
shape to the Kelvin- Helmholtz vortices [5 . Note that in 
the density cavities of Figure 3 the local magnetic field 
is approximately 45° out of the presented plane. 



LoAv density ribs Avithin the cavities 

The three dimensional structure of the cavities is stud- 
ied in detail in this section. The cavities appeared as 
thin electron current layers in two dimensional Particle- 
in-Cell simulations and in the contourplot on the plane 
z = Lz/2 in Figure 3. Figure 4 panel a, shows two isosur- 
face plots of the electron density at time Qcit = 14.5. The 



FIG. 2. Magnetic field lines at time Qdt = 12.1 in panel a. 
The color represents the magnitude of the magnetic field, nor- 
malized to the guide magnetic field Bq. Panels b and c show 
the reconnection electric field, Ez, normalized to BqVa/c from 
two different points of view at time Qcit = 12.1. 



electron density is normalized to the background density 
n5 = 0.1 no- The isosurface ng = 0.65 n5 (orange color) 
shows the cavity layer, approximately 1 di thick; the iso- 
surface ng = 0.4 n^ (grey color) reveals that regions with 
further lower density are embedded in the cavities. These 
structures are similar in shape to ribs, and for this rea- 
son we call these regions low density ribs. Separate low 
density ribs show nearly identical features. They have 
approximately equal size and bend along the magnetic 
field lines directions, close to the reconnection X line. 

Figure 4 panel b shows the density isosurfaces for de- 
creasing electron density values at successive times: the 
grey, red, blue and green isosurfaces represent regions 
with He = 0.45 nb at time Qcit = 12.1, ng = 0.3 n^ at 
time Qcit = 13.3, ng = 0.16 n^ at time Qcit = 14.5 
and ng = 0.15 n^ at time Qcit = 15.7 respectively. 
Two main results are found analyzing Figure 4 panel 
b. First, the cavities are perturbed by density waves 
that are almost stationary in space, far from the X line 
(12 di < X < 15 di) in a time period 12.1 < ftcit < 15.7. 
Closer to the X line (10 di < x < 12 di), the density waves 
start to overlap at different times, probably because of re- 
connection non-stationary dynamics and of the temporal 
variation of the electron velocity near the X line. Second, 
the density in the low density ribs decreases in time. Note 





FIG. 3. Contour plots of the electron density in the region x = 
4.8 di — 14.5 di,y = 4 di — 11.2 di,z = Lz/2 at different times. 
The plot shows the formation of the density cavity layers (blue 
regions), and the successive development of density ripples, 
reminiscent of Kelvin-Helmholtz vortices, along the cavities. 
The local magnetic field in the cavity is approximately 45° 
out of the presented plane. 



that the electron population can develop areas with den- 
sity as low as ne = 0.2 n5 when the ambient background 
plasma starts to dominate the separatrix edge. 

An intense parallel electron current develops along the 
cavities and contributes to the formation of the low den- 
sity ribs by progressively decreasing the density in lo- 
calized channels. This mechanism was identified previ- 
ously using two dimensional Particle-in-Cell simulations 
in Refs. |71|8], and our study confirms that electron ac- 
celeration inside the cavities is the main mechanism for 
creating the cavity density depletion. Figure 5 shows 
electron current streamlines in the cavity regions at time 
Qcit = 14.5. The electron current streamlines are orga- 
nized in bundles directed along the cavity channels. The 
electron velocity reaches approximately the peak 20 Vac, 
where Vac is the electron Alfven velocity. 

The low density ribs are aligned with magnetic field 
lines. The two dimensional studies performed in Ref. [5] 
reveal that the unstable electron Kelvin-Helmholtz waves 
grow rapidly in a plane perpendicular to the local mag- 
netic field. Figure 6 shows the magnetic field lines (red 
tubes): they are parallel to the grey Ue = 0.36 n^ iso- 
surface, which identifies the low density ribs, at time 
ftcit = 14.5. The low density rib structures are supported 
by the strong electric fields. An isosurface plot of the per- 
pendicular electric field intensity (E^ =E — E-B/|B|), 
shows the presence of a perpendicular electric field that 
oscillates in space in proximity of the cavities. The 
perpendicular electric field oscillates in space between 
^ —1.5 BqVa/c and ^ 1.0 BqVa/c values, with an os- 
cillation wavelength around one Larmor gyro-radius. Its 
maximum intensity is approximately 4.5 times the recon- 
nection electric field E^ (Figure 2 panels b and c). 




nne = 04nb 
□ ne = 0,65nb 



^^^^--..^^ 




-^m 




ffV 






Dcu 


= 12.1 


s^^ 


nn.= 


= 0.45 nb 


■ Oclt 


= 13.3 


f^mi 




^^ 


■ n.= 


= 0.3 nb 


■ Ocit 


= 14.5 


S^Il 


vw 


SW** 


■ n.= 


= 0.16nb 


□ Qcit 


= 15.7 


m 


1 


1 


an.= 


= 0.15nb 


^^ 




'■H..Mi 


l?^; 




\.^^^ 



FIG. 4. Panel a shows the two electron density isosurface 
plots for n = 0.65 nt (orange color) and n — 0.4 Uh (grey 
color) at time Q^cit = 14.5, revealing the cavities and the 
low density ribs within them. Panel b shows the electron 
density isosurface plots at successive times {Q^dt = 12.1 in 
grey, Q^cit = 13.3 in red, Q^cit = 14.5 in blue, and Q^dt = 15.7 
in green) for decreasing electron densities {ue — 0.45 Uh in 
grey, rie = 0.3 rih in red, rie = 0.16 rih in blue, and rie = 
0.12 Uh in green). The location of the low density ribs is 
approximately fixed in time, and their density progressively 
decreases. 



Bipolar electric field structures 

The carried-out simulations confirm the presence of 
bipolar parallel electric field {En = E • B/|B|) sig- 
natures along the cavities, reported by previous simu- 
lation studies J5l |29l [30]. Figure 8 shows an isosur- 
face plot oi Ell = ±0.3 BqVa/c with density isosur- 
faces He = 0.4 n5,0.65 n^ from two different points of 
view. Note that the net parallel electric field, that de- 
termines the electron acceleration inside the cavities, is 
much weaker than the bipolar electric fields and it does 
not appear in the plot. The presence of bipolar parallel 
electric field structures along the cavities, appearing as 
red and blue hemispheres, is evident. The plot reveals 
that the bipolar structures are isolated spherical forma- 
tions. The peak intensity of the bipolar parallel electric 
field signatures is of order 0.3 BqVa/c. The combination 
of electron density and parallel electric field plots in Fig- 
ure 8 shows that the bipolar electric field structures are 
elongated in the direction of the low density ribs, and of 
the magnetic field lines. Solitary electrostatic waves tend 
to appear between consecutive low density ribs in rela- 
tively higher density regions of the cavities. The positive 
hemisphere of the bipolar electric field precedes the neg- 
ative part in the direction toward the X point along the 
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FIG. 5. Electron current streamlines in the cavity region at 
time Qcit — 14.5. The color indicates the Jy value, normalized 
to eriQC. The electron current streamlines are organized in 
bundle in the cavity region. The X line is along the z direction 
dX X — 1^ di^y — 7.5 di (dashed black line). The streamlines 
at X < 10 di (left side in the x direction) and at x > 10 di 
(right side in the x direction) originate from seed lines atx = 
6.5 di^y = 8.5 di and x = 13.5 di^y = 6.5 di. The arrows 
indicate the electron flow direction in the x — z plane far from 
the X line. 




FIG. 6. Electron density isosurface plot for rie = 0.36 n^ with 
magnetic field lines in red tubes at time Q^dt — 14.5. The low 
density ribs are aligned with the magnetic field lines. 
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FIG. 7. Perpendicular electric field intensity isosurface plot at 
time Q^cit — 14.5. The perpendicular electric field reaches high 
intensity values (approximately 4.5 the reconnection electric 
field) and it oscillates in space supporting the low density rib 
structures. 
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separatrix. Panel b reveals that the the bipolar electric 
field structures form between the cavities and the outflow 
plasmas. 



DISCUSSION 

The reported Particle-in-Cell simulation shows the for- 
mation of the cavity layers and of low density ribs within 
them along the reconnection separatrices. The low den- 
sity ribs develop along the magnetic field lines, and their 
location remain approximately fixed in time. Their den- 
sity decreases progressively in time as electron currents 
evacuate them. A strong perpendicular electric field sup- 
ports the low density rib structures. 



FIG. 8. Electron density isosurface plot for n — 

0.4 rib, 0.65 Tih in orange and grey colors with parallel electric 
field isosurface plot for E/ / — ±0.3 BqVa/c in red and blue 
colors at time Q^dt — 14.5 from two different points of view. 
The bipolar parallel electric field structures develop along the 
low density ribs and magnetic field lines direction between 
the cavities and the outflow regions. Supporting material is 
provided online to analyze this configuration from different 
points of view. 



The origin of the three-dimensional low density ribs 
within the cavities is still unclear and further research 
is needed to identify the causes. Previous studies of 
guide field magnetic reconnection simulations [15 and 
a carried out spectral analysis of the wave in the cavity 
regions show that the frequency of the observed waves 



are in the range of the lower-hybrid frequency ujlh = 
cjpi/wl +Co'|g/r^^g, where Upe is the plasma frequency 

and (ice the electron gyro-frequency. Spacecraft observa- 
tions confirm the presence of lower hybrid waves also [8] . 
It has been suggested that lower-hybrid waves can be 
driven by electron beams [15] [3QH32] , or by density gradi- 
ents, as in the lower- hybrid drift instability [33]. The den- 
sity ripples in the cavities have been previously detected 
in two dimensional Particle-in-Cell simulations, and spec- 
tral analysis suggest that the instabiity is related to the 
electron MHD Kelvin-Helmholtz mode [5 . An additional 
possibility is that observed three-dimensional structures 
of higher density regions in proximity of the low density 
ribs are Alfven vortex filaments [34l [35] , driven by drift 
kinetic Alfven waves. The presence of Alfven vortices 
have been observed in cusp region [36] and in proximity 
of a magnetic reconnection X lines [371 135 • 

Moreover, the parallel electric field exhibits localized 
bipolar signatures in proximity of the cavity channels. 
These structures might be the signature of streaming in- 
stabilities such as the Buneman and/or two-stream in- 
stabilities [5], or originate from the beam-driven lower 
hybrid instability [30]. A further study is required to 
identify the origins of the bipolar parallel electric field 
structures. It has been found that bipolar electric field 
structures are approximately located in the cavities be- 
tween the low density ribs. One suggestive hypothesis 
is that the relative higher density regions in the cavities 
might act as waveguides, allowing the transport of bipo- 
lar electric field structures in channels between the low 
density ribs. 



CONCLUSIONS 

This paper presented a simulation study of three di- 
mensional structure of the cavity during collisionless 
magnetic reconnection with guide field. This study shows 
the presence of further lower density regions within the 
cavities. These low density ribs are directed along the 
magnetic field lines and they are supported by an in- 
tense perpendicular electric field. Bipolar parallel elec- 
tric field formations have been detected in the cavities 
between the low density ribs. Further research is nec- 
essary first to identify the instability or spectrum of in- 
stabilities leading to the low density rib structures and 
second to study their long-term evolution. In particu- 
lar, the use of open boundary conditions in the outflow 
direction is necessary to avoid the plasma recirculation 
at the boundary, enabled by the periodic boundary of 
the presented simulation. The open boundary condition 
will allow to follow the evolution of the low density ribs 
over a large period of time. The existence of the low 
density ribs within the cavities, of perpendicular electric 
fields, and of bipolar electric field structures is very im- 



portant the future NASA's Magnetospheric Multi Scale 
(MMS) mission. In particular, the reported features of 
guide field reconnection might be used to identify the 
proximity to the diffusion region. In fact, it has been 
suggested recently in Ref. [29 that bipolar electric field 
signatures can be used to detect the diffusion region and 
to control the spacecraft towards magnetic reconnection 
site. This research suggests that additional signatures, 
such as the low density ribs and associated perpendicu- 
lar electric fields, can possibly used as flag to detect the 
reconnection site also. 
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